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The Heyd-Scuseria-Ernzerhof (HSE) density functionals are popular for their ability to improve the accuracy 
of standard semilocal functionals such as Perdew-Burke-Ernzerhof (PBE), particularly for semiconductor 
band gaps. They also have a reduced computational cost compared to hybrid functionals, which results 
from the restriction of Fock exchange calculations to small inter-electron separations. These functionals are 
defined by an overall fraction of Fock exchange and a length scale for exchange screening. We systematically 
examine this two-parameter space to assess the performance of hybrid screened exchange (sX) functionals 
and to determine a balance between improving accuracy and reducing the screening length, which can further 
reduce computational costs. Three parameter choices emerge as useful: "sX-PBE" is an approximation 
to the sX-LDA screened exchange density functionals based on the local density approximation (LDA); 
"HSE12" minimizes the overall error over all tests performed; and "HSE12s" is a range-minimized functional 
that matches the overall accuracy of the existing HSE06 parameterization but reduces the Fock exchange 
length scale by half. Analysis of the error trends over parameter space produces useful guidance for future 
improvement of density functionals. 

PACS numbers: 31.15.eg,71.15.Mb 
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I. INTRODUCTION 

The exact, nonlocal form for the many-electron Fock 
exchange energy is known from Hartree-Fock theory. 
However, it is a standard practice in density functional 
theory (DFT) to compute this energy by integrating a 
local energy density per electron that is specified by the 
local electron density and its derivatives. One of the most 
popular of these semilocal density approximations is the 
Perdew-Burke-Ernzerhof (PBE) modet^. Accuracy can 
be improved by mixing the PBE exchange energy with 
a fraction of the exact nonlocal Fock exchange energy, 
producing hybrid functionals such as PBEO 2 . Exchange 
mixing can also depend on distance: the Heyd-Scuseria- 
Ernzerhof (HSE) density functional^ retains only short- 
range Fock exchange and preserves the accuracy of PBEO 
while avoiding the cost and pathologies^ of long-range 
Fock exchange. Unlike in traditional Kohn-Sham theory, 
hybrid functionals generate nonlocal potentials, which is 
described by a generalized Kohn-Sham (GKS) theory 5 -. 
Ideally, an exact GKS functional can produce both a 
highest occupied molecular orbital (HOMO) and a lowest 
unoccupied molecular orbital (LUMO) that have orbital 
energies equivalent to the negative ionization potential 
(IP) and electron affinity (EA), respectively. In practice, 
HSE and other short-range Fock exchange functionals 
are accurate for band gaps (IP+EA) of semiconductors^, 
but fail significantly for large-gap insulators 7 , molecules^, 
interfaces^, and nanostructuresi^. 

The HSE functional form defines a 2-dimensional space 



of DFT functionals-, set by the fraction of Fock exchange, 
a, at zero electron separation and a length scale, on 
which the short-range Fock exchange is computed, 



£PBE,LR (w)+ £PBE_ 



(1) 



The short-range Fock exchange (in hartrees) is calculated 
using the spinful Kohn-Sham density matrix p a a >(r, r'), 
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a ' Electronic mail: godotalgorithm@gmail.com 



while the long-range and remaining short-range exchange 
are derived from the exchange-hole formulation of PBE 11 . 
The intent of HSE was to achieve accuracy equivalent to 
PBEO, therefore the exchange fraction was limited to its 
PBEO value of a = 0.25. The HSE06 reparameterization 
of the HSE formi 2 - was also based on a variation of u> and 
not a. The remaining (uj, a) space has been explored only 
sparsely^ with results suggesting that different choices 
of uj and a improve the accuracy of different physical 
properties. 

HSE is not the first density functional to combine a 
short-range fragment of Fock exchange with a semilocal 
model of long-range exchange. The concept has appeared 
previously as screened exchange (sX) by Bylander and 
Kleinmani^, where a combination of short-range Fock 
exchange and the local density approximation (LDA) of 
correlation and long-range exchange, denoted 5 sX-LDA, 
was used to improve LDA band gaps of semiconductors. 
This application of short-range Fock exchange to band 
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FIG. 1. Approximation of the exponential exchange screening 
function with a scaled complementary error function. Both 
functions are displayed in the upper panel and their difference 
is displayed in the lower panel. 



the overall error over multiple test sets. The result is a 
modest improvement over the HSE06 parameterization, 
which comes from using a larger Fock exchange fraction 
(a = 0.313). Finally, we search for a range-minimized 
"HSE12s" parameterization that minimizes the value of 
uj^ 1 while preserving the accuracy of HSE06. This results 
in reduction of w _1 by a factor of two, which can reduce 
the cost of evaluating Eq. by reducing the number 
of p off-diagonals that need to be computed to achieve 
a target accuracy. After evaluating multiple tests sets, 
we observe significantly different behavior in two quanti- 
ties: molecular IPs and EAs as approximated by HOMO 
and LUMO energies. We can explain this phenomenon 
as a fundamental limitation of the HSE functional form 
that results from the use of an environment-independent 
and homogeneous screening of exchange that is unable to 
account for inhomogeneous screening environments. 



II. COMPUTATIONAL DETAILS 



gap estimation is motivated by the Coulomb-hole-plus- 
screened-exchange (COHSEX) approximation to electron 
quasiparticle theory 1 - in combination with the Thomas- 
Fermi model of screening. Based on this motivation, an 
exp(— r/rxF) screening factor is used in place of erfc(cjr) 
in Eq. ^ and the Fock exchange fraction is set to unity 
(a = 1). The difference in screening form is insignificant, 
since the two variants can approximate each other, 

exp(-a;) « 0.95erfc(0.58x), (3) 

with a maximum pointwise error of ps 0.05 (see Fig. [TJ). 
However, HSE uses a single uj parameter for all systems, 
while the t"tf parameter in sX-LDA is set to the Thomas- 
Fermi screening length based on average valence electron 
density ~p (in atomic units), 



System-dependence of the Fock exchange length scale is 
the most important difference between HSE and sX-LDA, 
especially because p is not well-defined for a molecule 
without an arbitrary definition of molecular volume. 

Given the flexibility of the HSE form, we fully explore 
the (uj,a) parameter space to determine if its accuracy 
and efficiency can be improved further. Also, we assess 
the variation in results over this space to explain why 
HSE is accurate and reveal its unresolved shortcomings. 
First, we search for a "sX-PBE" functional in HSE space 
that reproduces the performance of sX-LDA, which we 
define by setting a — 0.95 and optimizing uj for band gap 
estimation. Despite the use of a system-independent u) 
instead of a system-dependent ttf, sX-PBE is able to 
surpass the accuracy of sX-LDA. This can be explained 
by the small degree of variation in ttf values. Next, we 
search for a "HSE12" parameterization that minimizes 



The analysis we present in this paper is based on DFT 
calculations of isolated molecules and periodic solids. All 
molecular calculations are performed in GAUSSIAN 09^. 
All periodic solid calculations are performed in a modified 
version of VASP 5.2il~— . All final computational results 
not directly reported in the paper, as well as all molecular 
and crystal structures used in this study, are compiled in 
the supplementary material^!. 

All GAUSSIAN calculations are performed using the 
6-311++G(3df, 3pd) basis set, consistent with previous 
HSE studies^. The exchange hole form of PBE±i that 
is used within HSE is applied to all calculations, even 
when w = 0. Default settings are used except for the 
fraction of non-convergent calculations, which we fix by 
switching to the quadratically-convergent self-consistent 
field (SCF) method. 

All VASP calculations are performed using the manual- 
recommended PBE projector augmented wave (PAW) 
pseudopotentials-- with their default planewave basis 
sets, which are assumed to be transferrable to all DFT 
functionals considered in this paper. The Brillouin zone 
integration is performed with the tetrahedron method 
on a T-centered grid, 12 x 12 x 12 for cubic solids and 
12 x 12 x 8 for hexagonal solids. The Fock exchange is 
down-sampled by half in each direction, consistent with 
convergence studies—. Fully sampled Fock exchange is 
applied as an additional perturbative correction to the 
band gapS4. Minor modifications to the current version 
of VASP are required to calculate the 3 sX-LDA variants 
considered in this paper—. 

Comparing DFT results to experiments may require 
finite temperature and quantum nuclear corrections that 
are not typically included within DFT itself. The largest 
corrections occur for formation energies. We account for 
these corrections at the level of G3 theory for the G3/99 
test set>2£, but omit them for all other test sets. As DFT 
becomes more accurate, the importance of accounting for 
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these effects to correctly assess accuracy will grow. Also, 
spin-orbit coupling is omitted from all calculations. This 
can have a large effect on band gaps, but the SC/40 test 
set£ that we utilize has removed large spin-orbit effects 
from its experimental values. 

Error distributions are quantified using the mean error 
(ME) with a (theory - experiment) sign convention, the 
mean absolute error (MAE), and the root-mean-square 
error (RMSE). Figures show MAE, which is at present 
the preferred error average in quantum chemistry. 

III. RESULTS AND DISCUSSION 

The physically relevant region of the HSE parameter 
space is < to < oo and < a < 1. For sampling and 
visualization purposes, the uj parameter is mapped to a 
bounded and dimensionless range separation parameter, 

~ 2 / u \ 

uj = — arctan , (5) 

7T \^HSE06/ 

with whseo6 = 0.208 A -1 and < uj < 1. This aligns 
the HSE06 parameterization to the center of the range 
separation axis, 5hseo6 = 0.5. All points in HSE space 
are specified as {uj, a) ordered pairs. In some cases, the 
range separation parameter is given as a screening length, 
which we derive using Eq. (£[]) as ttf = 0.58/w. 

A mixture of molecular and bulk solid tests are used to 
assess the accuracy of HSE functionals. The SC/40 sel£ 
contains 33 band gaps and 42 lattice constants of binary 
and elemental semiconductors (lattice constants of the 3 
hexagonal structures are omitted from our study). The 
G3/99 se1j2& contains 223 formation energies, 86 IPs, and 
58 EAs for small molecules (excited states of SH^ and 
are omitted from our study). The BH42/04 set22 
contains 39 distinct barrier heights mostly for pairwisc 
interactions between small molecules. The T-96R set2& 
contains 96 bond lengths of small molecules. All of these 
test sets have been used before to assess the accuracy of 
HSE functionals^^. 

A. SC/40 band gaps 

The comparison between sX-LDA and HSE is confined 
to the SC/40 band gaps because most sX-LDA studies 
have focused on band structures. We connect sX-LDA to 
HSE through a sequence of four functional modifications 
with the results shown in Table |U In sX-LDA, the long- 
range exchange energy is computed from integration of a 
globally- weighted local exchange energy density (LDAq), 

£ x LDAg ' LR = / dvp{r)e^[p{v)]f (pV3 rTp ) . (6) 

The use of p instead of p(r) in the weight function / was 
initially chosen based on improved band gaps of silicon^. 
/ is based on an exchange-hole formulation of LDA, and 



TABLE I. Comparison of ttf (from Eq. (Q, without semicore 
d-states in p) and band gaps of the SC/40 sel£ using 5 sX-DFT 
functionals described in the text, labelled by their exchange- 
correlation model, screening function, and screening length. 
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consistency with the exchange-hole formulation of PBE 
used in HSE 11 requires a local weighting (LDAl), 

£ x LDAl ' LR = / dvp{v)e^[p{v)]f (p(ry/ 3 r TF ) . (7) 

The effect of switching between LDAq and LDAl in sX- 
LDA is shown in the first two sX-DFT columns of TableQ] 
As in silicon, this change results in a systematic band gap 
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33 SC/40 band gaps (MAE) 
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FIG. 2. Semiconductor band gap errors over HSE space. A 
subset of functionals with approximately equivalent accuracy 
extend from (0,0.19) to (0.83,1). 



underestimation in relation to the more accurate LDAq ■ 
The following two modifications are replacement of the 
exp(— r/rxF) screening function with an optimal erfc(wr) 
approximant as in Eq. ^ and the replacement of the 
LDA model of correlation and long-range exchange with 
the PBE model used by HSE. As shown within the third 
and fourth sX-DFT columns of Table HI this produces 
small changes in the band gaps. The final modification 
is to choose a single uj value for all materials to recover 
the HSE form. The variations in ttf appear to be small 
because of its sixth-root dependence on 1/p in Eq. £1|). 
As a result, a single optimized uj value is able to surpass 
the performance of sX-LDA with material-dependent ttf 
values for this test set. We denote the optimized HSE 
approximation to sX-LDA as sX-PBE. 

The consistency of our results with previous sX-LDA 
results is mixed. We agree with prior result o 5 ' 14 ' 29 ' 30 on 
Si but produce significantly smaller band gaps for Ge and 
GaAs compared to other reported sX-LDA values 5 -^ - — . 
These discrepancies can be explained by the absence of 
valence-core interaction terms^ in prior sX-LDA results. 
The changes are small in Si but large in Ge and GaAs. 
vasp includes the effect of these missing terms as all- 
electron corrections within atomic spheres^. 

The MAE of SC/40 band gaps over the whole (uj, a) 
HSE space is plotted in Fig. Errors near (1,0) come 
from band gap underestimation in the PBE model that 
is largely attributed^ to the absence of discontinuities 
in SE^ E /Sp(r). Errors near (0, 1) come from band gap 
overestimation in Hartree-Fock theory^ (the additional 
correlation potential has relatively little effect here). In 
a GKS theory context, the 5 eV error variation in going 
from (1,0) to (0,1) originates from tuning the amount of 
nonlocality in the GKS potential to compensate for the 



lack of 8E xc /5p(r) discontinuities in the semilocal model. 
Within the HSE space, there is a 1-dimensional subset 
of similarly good functionals for band gap estimation as 
judged by MAE. This non- uniqueness is not surprising 
given that GKS theory introduces greater flexibility in 
functional form without adding a commensurate amount 
of new mathematical or physical constraints. Requiring 
that GKS theory reproduce the exact 1-electron reduced 
density matrix (1RDM) instead of just electron density 
can achieve uniqueness^. However, this introduces the 
ongoing problem of constructing GKS functionals that 
produce the fractional occupation numbers characteristic 
of lRDMs for interacting electrons, even for pure states 
with an integer number of electrons^. 

We can argue for the accuracy of band gaps in HSE 
based on its similarity to an electron quasiparticle theory. 
Specifically, the simplification of Hedin's equations using 
the COHSEX approximation 15 results in the self-energy 
(assuming a natural spin-dependent extension) 



E coHs EX(ry) = _^ CT , (r y )w ^, (r y) 



(8) 



with a non-local sX operator and a local potential derived 
from the Coulomb hole, both determined by a statically 
screened Coulomb interaction W a , a i (r, r'). Minimal, yet 
realistic, semiconductor screening models^ contain the 
bulk static dielectric constant £rj and a Thomas-Fermi- 
like screening length ttf- A simplified example is 



HW(r,r') 



(l-e - 1 )exp(-|r-r'|/r TF ) 



(9) 



which has a constant Coulomb hole, (e^ 1 — l)/(2rxF), 
corresponding to the polarization energy of an electron 
in a spherical cavity of radius ttf within a dielectric of 
permittivity e^r^. We presume that oversimplification of 
W breaks the connection between sX and the Coulomb 
hole, the latter being reasonably approximated by the 
PBE local potential that corresponds to correlation and 
long-range exchange. 

The position and shape of the 1-dimensional subset of 
accurate functionals within Fig. ^fy can be explained by 
the connection between HSE and COHSEX. For any one 
semiconductor, we have one constraint, the band gap, 
and two free parameters, (uj, a), therefore each can be fit 
perfectly along a line in HSE space. The near-overlap 
of all these lines results from similar parameters for all 
members of the SC/40 set. The effective W of HSE, 



W™?(r,r>) 



erfc(u;|r — r'|) 



(10) 



can fit Eq. © when r^p = or e = 0. We assume 
similar gap-invariant variations of (r^^E^ 1 ) in Eq. Q 
from reasonable physical values to (r^p = 0,e^ > Eq 1 ) 
or (rv^ F < t^. f ,Eq = 0). This explains the optimal a sa 
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FIG. 3. Molecular formation energy errors over HSE space. A 
subset of functionals with approximately equivalent accuracy 
extend from (0,0.3) to (0.93,1). 



Eq « 0.2 that is larger than the average value for the 
SC/40 set, £q 1 ss 0.1. It also explains the behavior in 
Table Q] of an optimal sX-PBE screening length that is 
longer than the Thomas-Fermi screening lengths. This 
malleability also explains how the same arguments can be 
used to support both hybrid- and sX^ DFT functionals, 
with comparable validity and empirical success. 



B. G3/99 formation energies 

The MAE of G3/99 formation energies over HSE space 
is shown in Fig. ([3]). As with the SC/40 band gaps, there 
is a 1-parameter family of similarly accurate functionals, 
here with a larger exchange fraction (30% versus 20%) in 
the hybrid limit and a smaller effective screening length 
(0.33 A versus 0.79 A) in the sX limit. A negligible in- 
dependence for screening lengths larger than 2.4 A (S < 
0.55) shows that the PBE model of long-range exchange 
performs well beyond that length scale, at least within 
small molecules. GKS orbitals are only weakly dependent 
on details of the exchange functiona l 40 ' 41 , therefore errors 
in the exchange functional itself are the main source of 
MAE variation. In particular, the largest errors occur 
near (0, 1), where the exact exchange energy is used. This 
is an example of error cancellation between exchange and 
correlation models that is often observed in DFT—. 

The errors in formation energy can be attributed to a 
systematic bias. Similar to the results for band gaps, the 
formation energy is underestimated near (1, 0) in Fig. ^ 
and overestimated near (0, 1), with ME passing through 
zero near the region of minimal MAE. Here, we use the 
sign convention of positive formation energies for stable 
molecules. We observe that the PBE exchange energy 



overestimates the Fock exchange energy, > i 

with larger errors for molecules than for atoms. Also, we 
observe that the PBE correlation energy underestimates 
the experimental correlation energy, E^ BE ^ E° xp , with 
errors again larger for molecules than for atoms. These 
observations are consistent with a cancellation of errors 
between exchange and correlation contributions to the 
formation energy. 

The similar error trends in Figs. © and (J3J suggest 
that HSE functionals are "right for the right reasons". 
The ability of HSE functionals to estimate band gaps is 
explained by quasiparticle theory, which can be extended 
to cover total energy approximation. This is related to 
ongoing efforts^ to compute accurate total energies from 
self-consistent 1-electron Green's functions. Specifically, 
the HSE functional can be interpreted as the Galitskii- 
Migdal formula 44 applied to the frequency- independent 
COHSEX self-energy in Eq. © to calculate an exchange- 
correlation energy, 

E ^ = \Y.J drdrV, CT (r',r)£ ff , CT <(r,r'). (11) 

(7, a 1 

The effective HSE self-energy is particularly suited for 
total energy calculations since it is constrained to exactly 
reproduce the total energy of the uniform electron gas 
and PBE gradient corrections thereof. 

C. G3/99 ionization potentials and electron affinities 

Vertical IPs and EAs can be computed with a GKS 
functional using either total energy differences between 
charge states (ASCF) or the HOMO/LUMO eigenvalues, 
if the functional is accurate and free of discontinuities in 
8E xc /5p(r). The G3/99 test set contains adiabatic IPs 
and EAs, which requires the use of relaxed geometries 
in ASCF calculations and a relaxation energy correction 
for eigenvalue-based vertical excitation energies (the total 
energy difference of the neutral molecule between neutral 
and ionized geometries). ASCF IPs and EAs are shown 
in Figs. (|4ji) and (0Ji). The variation of errors is smaller 
than in the formation energies of Fig. ([3]), which results 
from increased cancellation of errors because two charge 
states of a molecule are more similar than a molecule and 
its dissociated atoms. There is still the same subset-of- 
accurate-functionals trend, but it does not stand out as 
strongly without a large background of systematic errors. 
The eigenvalue-based IP and EA estimates in Figs. (HJd) 
and ((4J3) show much greater variation over HSE space, 
comparable to the band gaps in Fig. ©. There is still 
a subset of accurate functionals, but it is shifted to large 
Fock exchange fractions (a > 70%) that are incompatible 
with accurate formation energies or band gaps. 

The discrepancy between ASCF and eigenvalue-based 
IPs and EAs constitutes both delocalization error— and 
size-consistency error—. Compared to ASCF values, the 
eigenvalues underestimate IP and overestimate EA near 
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FIG. 4. IP and EA errors in HSE space, comparing ASCF to eigenvalue-based estimates and including polarization corrections. 



(1, 0) in HSE space and vice versa near (0, 1). If a hole 
or electron is added to a dilute gas of identical molecules, 
the charge will delocalize to avoid spurious Coulomb self- 
interaction and approach the erroneous eigenvalue-based 
IP/EA for functionals near (1,0). The opposite effect, 
localization of charge onto a single molecule, occurs near 
(0, 1) and forces IP/EA to remain at the more accurate 
ASCF value. This dichotomy is demonstrated in Fig. ([5]) 
for a dilute He ra gas with n = 1, • ■ • , 20. DFT functionals 
with derealization behavior are not size consistent, and 
making use of the ASCF accuracy requires a single well- 
defined charge center such as a small molecule or a point 
defect inside a crystal^. The subset of functionals that 
are accurate for formation energies and band gaps is well 
within the derealization regime and uniformly lacks size 
consistency. In the localization regime, the eigenvalues 
are no longer equivalent to charge excitations because 
the localized charges produces nonlinear corrections. 

Considering IP as an example, we can account for the 
difference between ASCF and eigenvalue-based estimates 
with a continued appeal to quasiparticle theory, building 
on the arguments for band gap accuracy in section 1111 Al 
and total energy accuracy in section UlIBI For a system 
X and a DFT functional with a COHSEX-like form, the 
IP discrepancy can be written as 

E(X+) - E(X) = -enoMO - A relax 
+ J drdr'(^-l^-^(r,r'))p+(rV+(r') 

(12) 



by separating orbital relaxation effects A re i a x from terms 
that assume the cation to have the same W, orbitals, 
and energies as the neutral system. Within Hartree-Fock 
theory, the last term vanishes and this is the essential 
content of Koopmans' theorem 4 ^, which guarantees that 
ASCF doesn't produce a larger IP than — ehomo since 
Areiax > 0. The last term is a polarization correction 
to the self-interaction of the hole charge p + and, in an 
accurate quasiparticle theory, it should cancel A rc i ax to 
produce consistency between IP estimates. If we assume 
this to be true for a COHSEX model with a realistic W, 
then errors in HSE can be attributed to the lack of an 
unscreened Coulomb tail, W l/|r — r'|, extending out 
from a small molecule. The net effect is over- screening 
in small molecules compared to semiconductors, which is 
consistent with the offset between Figs. © and (0b). 

The difference between 1/ |r — r'| and W must originate 
from an electronically polarizable substance. Since VFhse 
remains screened even in vacuum, the effect there can be 
assigned to unphysical vacuum polarization. This error 
cannot be removed without environment-dependence and 
inhomogeneity in Whse, but we can model it by adding 
comparable errors into ASCF calculations. Specifically, 
we introduce an artificial vacuum polarization into total 
energies with the polarizable continuum model (PCM) 49 . 
Over-screening errors of eigenvalue-based IPs and EAs 
are corrected by subtracting out the errors modeled by 
the difference between ASCF+PCM and ASCF values. 
The results in Figs. (j4j:) and (j4f) are less accurate than 
ASCF values, but they succeed in recovering previous 
error trends in HSE space. Similar over-screening errors 
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FIG. 5. IP of a dilute He„ cluster, comparing — £HOMo(He„) 
to £(He+)-£(He„) using the HSE06 functional and full Fock 
exchange with PBE correlation (HF+c), which correspond to 
the points (0.5,0.25) and (0.0, 1.0) in HSE space. In HSE06, 
the ASCF and eigenvalue-based results converge for n — > oo. 
The electron-hole of He^ is localized on one atom with HF+c 
and uniformly divided among all n atoms with HSE06. Any 
n-dependence is a size-consistency error. 



should also occur for localized defect states in crystals, 
but with the discrepancy in polarization reduced by £q . 
Another solution is to include long-range Fock exchange 
to model the correct asymptotic of W, with its range 
parameters tuned to model the size of the molecule and 
approximately satisfy Koopmans' theorem^. Given our 
results with HSE, caution must be exercised when tuning 
range-separated hybrid functionals into a regime where 
HOMO and ASCF values of the IP match, because this 
may degrade the accuracy of other physical properties. 



D. Conformational tests: BH42/04 barrier heights, T-96R 
bond lengths, and SC/40 lattice constants 

The final set of tests are based on atomic conformations 
that are less pronounced than full formation, which leads 
to less variation of error in HSE space than the formation 
energies in Fig. ([3]) . For the barrier heights in Fig. (|6ji) , 
the optimal uj value at a = 1 is consistent with previous 
results The error trend in molecular bond lengths in 
Fig. (JSJd) matches well with molecular formation energies 
in Fig. (|3]). The increased w-dependence at small w and 
large a of the lattice constants in Fig. I^jp) are the result 
of unconverged Brillouin zone sampling of the long-range 
exchange tail. In all cases, the minimum in MAE closely 
corresponds to the zero of ME. The net effect is to reduce 
small systematic errors in PBE that overestimate bond 
lengths and lattice constants and underestimate barrier 
heights. 



E. Outlook on method development 

Our combined computational survey and analysis of 
HSE functional space suggests several new directions for 
method development with varying degrees of expected 
difficulty and efficacy. We have shown the approximate 
correspondence of sX-LDA to the HSE functional form. 
Given the limited performance of sX-LDA for properties 
other than band gaps^, we conclude that HSE is both 
conceptually and empirically superior because of its use 
of PBE instead of LDA for correlation and long-range 
exchange. Any benefit to the system-dependence of the 
sX-LDA screening length that is negligible for solids and 
ill-defined for molecules is offset by optimization of a 
system-independent HSE screening length. Within the 
HSE form, we can recommend new parameters that are 
modest improvements over HSE06. We base this on an 
aggregate error metric shown in Fig. ([7]), 



1 



tests 



# of tests 



MAEj(D,a) 



■ MAE,- 



- 1 



(13) 



The eigenvalue-based tests are excluded because of their 
systematic bias. This metric puts emphasis on the more 
sensitive tests and is zero if the MAE of each test can 
be simultaneously minimized. The "HSE12" functional, 
with uj = 0.185 A -1 and a = 0.313, optimizes accuracy 
by minimizing the metric. The "HSE12s" functional, 
with ui = 0.408 A -1 and a = 0.425, minimizes screening 
length while preserving the accuracy of HSE06. HSE12s 
has half the screening length of HSE06, which can reduce 
computational costs in implementations that make use of 
spatial decay in short-range Fock exchange. A summary 
of functionals and their performance is given in Table |H] 

The reduced cost of HSE12s can be demonstrated by 
the relative speed-up of an example: bulk aluminum. In 
GAUSSIAN (using the basis from Ref. 51), the relative 
speed of HSE12s compared to HSE06 is 1.6. The cause 
of this improvement is exchange integral screening 5 -^. In 
VASP, the speed-up is 1.4. Here, the improvement comes 
from reduced Brillouin zone sampling of the sX term for 
HSE12s compared to HSE06 ?3 >, requiring 5x5x5 points 
instead of 6 x 6 x 6 to achieve equivalent convergence. 

Further progress in increasing functional accuracy will 
require alterations of the HSE form. More complicated 
screening forms have been proposed 5 ^, which enable more 
realistic W models, such as in Eq. ©, and may help to 
align the optimal choice of screening for multiple physical 
properties. The PBE model was designed for exclusively 
semilocal treatment of exchange, and a recent result 53 
shows the benefits of refitting PBE dependent on how 
exchange is split into Fock and semilocal components, 
because it enables a fit of PBE to the residual, rather 
than total, exchange energy. Based on the strong ties 
between HSE and COHSEX outlined in this paper, it 
may be worthwhile to pursue semilocal DFT functionals 
based on a sX/Coulomb-hole partition^! instead of an 
exchange/correlation partition. Given the observed error 
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FIG. 6. Binary reaction barrier height, molecular bond length, and semiconductor lattice constant errors in HSE space. 



trends in HSE space, it may be possible to construct a 
scmilocal sX functional by taking the subset of accurate 
functionals to the (u, a) — > 00 limit. Unfortunately, these 
suggestions are not likely to fix the large errors in HSE 
HOMO/LUMO eigenvalues for molecules, because they 
do not introduce any nonlocal screening effects. 

The inevitability of nonlocality can be argued in the 
context of GKS theory as an approximation to COHSEX 
theory. The screened exchange operator in Eq. ([8]) is a 
direct product of p and W, which is a spatially localized 
operator for basic electronic phases. In an insulator 55 , 
p decays exponentially with increasing |r — r'|, while W 
decays only algebraically. In a metal, it is p that decays 
algebraically and W that decays exponentially. In either 
case, the product is localized by one constituent while the 
other builds nonlocal dependence into the quantitative 
details. Hypothetical examples of nonlocal dependence 
are electron interference in p resulting from proximity to 
a scatterer in a metal or increased screening in W from 
proximity to a highly polarizable object in an insulator. 
In the context of DFT, the exchange-correlation hole is 
observed to be spatially localized, but no limits have been 
placed on its sensitivity to distant perturbations 56 . 

DFT modeling can be roughly split between quantities 
that are tractable (albeit expensive) to compute exactly 
(e.g. Fock exchange energy) and those that cannot be 
generally computed (e.g. exact correlation energy). At 
the level of COHSEX theory, W is contained within the 
first category and can be explicitly constructed from an 
independent-electron response function, 



€i - £j 



W-ArA) = 



1 -1 



X° <(r,r'), 



(14) 



for orbitals ipi, energies q, and occupations and where 
1-1 ' refers to operator inversion. Using such forms for W, 
promising results for total energies and charge excitations 
have been demonstrated for self-consistent quasiparticle 
methods 5 -^— and much of this accuracy is preserved by 



the COHSEX approximation^ -. It is conceivable that the 
precision and reliability necessary to attain the goal of 
"chemical accuracy" will only be achieved with a detailed 
treatment of W in the same way orbital-free DFT has 
not been able to achieve a level of accuracy comparable 
to methods containing the details of electronic orbitals 61 . 
The construction of p from electronic orbitals, 



(15) 



contains its essential nonlocal character and accounts for 
much of the success of Kohn-Sham DFT. In principle, 
p and W have equal importance in COHSEX theory, 
and putting this into practice in COHSEX-inspired DFT 
functionals will require the use of Eq. (fl"4"f . 

The form of W in Eq. (fl~4")) highlights an open problem 
in self-consistent quasiparticle theory. Quasiparticle and 
GKS methods are designed to model charge excitations, 
but the electronic polarizability represented in Eq. (| 14[) 
is supposed to arise from neutral excitations. The use of 
quasiparticle energies in Eq. (fT4"|) produces a systematic 
overestimate of band gaps, which has been removed with 
vertex corrections encoding electron-hole interactions 6 ?. 
Otherwise, the non-self-consistent evaluation of Eq. ([14"]) 
using DFT methods is found to produce accurate mean- 
field models of static polarization 63 . Similar observations 
have been made for HSE: that it approximates the optical 
gap rather than the charge gap when there is a significant 
difference between the two 6 ^, even though it is a GKS 
theory that should be approximating the charge gap. A 
heuristic explanation for this effect is that the universal 
lack of exchange screening in HSE can be attributed to 
electrons (holes) being screened by both a ground-state 
polarization and an excess hole (electron). Just as with 
the excess polarization argument in section IIII C[ this 
energy can be removed from the eigenvalues to improve 
their correspondence to quasiparticles, and in this case it 

as the Coulomb 



can be re-identified for Clumo 



C HSE 
"HOMO 



binding energy between an electron and hole. It may be 
possible for a self-consistent COHSEX-based theory to 
avoid vertex corrections if it can define a consistent pair 
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FIG. 7. Aggregate error in HSE space. The minimum of this 
error metric occurs at the HSE12 parameterization. HSE06 
and HSE12s lie on a contour of equal error. 



of mean fields: a polarization mean field, meant to model 
neutral excitations for Eq. (|14p , and a quasiparticle mean 
field containing the COHSEX self-energy from Eq. (JHJ. 
Based on existing theory, it is possible to construct GKS 
and optimized effective potential (OEP) 65 methods from 
the same COHSEX-based total energy functional, with 
the GKS form serving as the quasiparticle theory and the 
OEP form serving as the polarization model. However, 
a consistent theoretical framework for such a "double" 
mean-field theory does not yet exist. The theory would 
need to specify the interrelations between the orbitals 
and energies of the two mean fields and how they should 
be used in concert to calculate a total electronic energy. 



IV. CONCLUSIONS 

The novelty of our study comes from HSE's precarious 
perch between empirical and non-empirical DFT, with 
just two tunable parameters. This enables a thorough 
sampling of parameter space, similar to recent studies 
on the PBE functional 66 . Besides the immediate benefit 
of further fine-tuning accuracy, the study reveals error 
trends supporting a connection between the HSE form 
and quasiparticle theory, which explains its success in 
approximating semiconductor band gaps. The degree to 
which multiple physical properties can be simultaneously 
optimized by a common set of parameters is pleasantly 
surprising. Unpleasant and unsurprising^ are the large 
discrepancies between ASCF and eigenvalue estimates of 
small molecule IPs and EAs, which we argue to originate 
from the absence of environment-dependent screening in 
the HSE functional resulting from discrepancies between 
W mE in Eq. (ID) and W in Eq. flU). It is encouraging 



TABLE II. Summary of important functionals in HSE space 
and their performance on various tests discussed in this paper. 
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that this error can be fixed by correcting the long-range 
tail of W— . However, this simple fix is not transferrable 
to more complicated systems such as a molecule near a 
metal surface, where a model W such as Eq. © is not 
flexible enough to set e^ 1 = for the metal half-space 
and Eq 1 = 1 for the empty half-space. There is much to 
be gained from incorporating more sophisticated models 
of W into the screened Fock exchange component of GKS 
functionals, such as Eq. (|14p or a reasonable facsimile 
thereof^. 
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